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SUMMARY 

A theoretical and experimental investigation has been 
made of the behavior of a cantilever beam in transverse 
motion when its root is suddenly brought to rest. Equations 
are given for determining the stresses, the deflections, and 
the accelerations that arise in the beam as a result of the 
impact. The theoretical equations, which have been confirmed 
experimentally, reveal that, at a given percentage of the 
distance from root to tip, the bending stresses for a 
particular mode are independent of the length of the beam 
whereas the shear stresses vary inversely with the length. 

INTRODUCTION 

When an airplane lands, the vertical component of the 
velocity is rapidly reduced to zero. In the absence of a 
thorough analysis of the stresses that arise from such 
shocks, it is customary for engineers to assume that the 
landing loads are static and independent of the elastic 
properties of the structure. As an initial step in the study 
of elastic structures under shock loads, an investigation has 
been made to determine the effect on a simple structure of 
the sudden arrest of its motion and the effect of the 
geometry of the structure on the stresses that result. The 
particular case treated in this report covers the basic 
problem of the instantaneous arrest of the root of a moving 
cantilever beam. The solution of this problem gives the 
energy consumed in exciting the different modes of vibration 
and the stresses, deflections, and accelerations that result 
throughout the beam. 

This investigation is based on the usual engineering 
beam theory in which the deflections are considered to be the 
result of bending alone and shear deflections are neglected. 
The theory, as applied to ordinary beams, gives reasonably 
good results as long as the distance between inflection 
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points is greater than a few times the depth of the beam. 
When this theory for beam action is used in vibration 
problems, such as the problem in the present paper, the 
results are satisfactory for those modes of vibration for 
which the nodes are not too close together. This report 
summarizes the results of a theoretical solution, given in 
the appendix, and presents an experimental verification 
of these results. 



SYMBOLS 

E modulus of elasticity 

Y weight density of material 

\ coefficient of equivalent viscous damping of 

material 

c velocity of sound in material 

g acceleration of gravity 

L length of beam 

I moment of inertia of cross section of beam about 

neutral axis 

A cross-sectional area of beam 

p radius of gyration of cross section of beam 

x coordinate along beam measured from root 

y distance from neutral axis of beam to any fiber 

t time, zero at impact 

p operator 



n integers 1, 2, 3, etc. designating a particular 

mode of vibration 

9 n nth positive root of 1 + cos 0 cosh 8 = 0 

:jo n undamped natural angular frequency of nth mode, 

radians per second Ipc - 
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V 

w(x,t) 
w n (x, t) 

a(x,t) 
a n (x,t ) 

a ( x , y , t ) 
T(x,t) 

T n (x,t) 

A n ,• * i 

B n 

Cn 



damped natural angular frequency of nth mode, 



When - — — > 1, the "frequency" Is defined 



velocity of beam prior to impact 

deflection of beam at station x and time t 

deflection of beam at station x and time t 
for nth mode of vibration 

acceleration of beam at station x and time t 

acceleration of beam at station x and time t 
for nth mode of vibration 

bending stress in beam at station x, distance 
from neutral axis y, and time t 

bending stress in beam at station x, distance 
from neutral axis y, and time t for 
nth mode of vibration 

average shear stress over cross section of beam 
at station x and time t 

average shear stress over cross section of beam 
at station x and time t for nth mode of 
vibration 

bending- stress coef f i cient 

shear- stress coefficient 

deflection coefficient 



radians r)er second 
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RESULTS AND CONCLUSIONS 
Theoretical 



When a cantilever beam under uniform translation in 
a direction perpendicular to its length has its root 
instantaneously brought to rest, there is excited a 
theoretically infinite number of modes of vibration. 
With each successive mode, damping has an increasing 
influence upon the frequencies and amplitudes of vibration 
and, for sufficiently high modes, even changes the type 
of motion from oscillatory to nonoscillatory motion. In 
the lower modes, however, damping has little effect, and 
only terms of the first order in damping need to be included 
in the equations. Only the equations applicable to the 
lower modes, which alone are of importance in any practical 
case, are presented in this section of the paper. For a 
more complete treatment of damping, see the appendix. 

The angular frequencies ( 2tt times the frequencies 
in cps) are given by the equation 

9 2 

^n = Pc -75T (1) 



where 9 n has the following values for successive modes 
of vibration: 

B 1 = 1.87510)4- 85 - lit. 137168 

6 2 = ^.69^098 9 6 as 17.278759 

65 = 7.85ll757 9 n * \ (2n-l)TT, n > 6 

8^ = 10.9955U 

The energy that the beam possesses before impact is 

consumed in exciting the various modes of vibration and 

is distributed among the modes as follows: 
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Mnrlp T» 
1V1 vj w w , 11 


Percentage of energy 


1 


61.5 


2 


18.8 




7.k 


j 


5-2 


1 


1.9 




1.3 


7 tO co 


6.1 



This distribution of energy among the different modes of 
vibration is presented graphically in figure 1. 

All stresses, deflections, and accelerations are 
damped sinusoidal functions of time and vary along the 
length of the beam. The bending stress a n (x,y,t) and 
the average shear stress T n (x,t), associated with the 
nth ?node of vibration, are given by the equations 

o n (x,y,t) = A n J £se " sin co n t (2) 

r n (x,t) = B n g ^Ee~ sin fj0n t (3) 

The variation of the dimensionless coefficients A n 
and B n with x/L is given for n = 1, 2, and 3 in 
figures 2 and 3. The highest values of A n and 
and hence the highest stresses, occur at the root of the 
beam. These values, for the first six modes, are 



Mode, n 


A n at root 


E n at root 


1 


I.566 


2.11+6 


2 


.863 






-509 


3.99k 


I 




k.oo 




.283 


i+.oo 


I 


.231 


k.OO 



The foregoing values of A n and Bn at the root are 
presented graphically in figure 

The maximum values with respect to time of a n (x,y, 
and "T n (x,t) associated with the nth mode of vibration, 
when the effects of damping are neglected, are 
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a n (x,y) = A 




(k) 



T n(x) = 3 n | £e 



(5) 



The deflections w n (x,t) for the nth mode of 
vibration are given by the equation 




p n 



(6) 



The accelerations a n (x,t) for the nth mode, when 
damping is sufficiently small, are given by 



a n (x,t) = -GO 



n 



w n (x,t ) 



(7) 



The variation of the dimensionle ss coefficient C n 
with x/L is given for n = 1, 2, and J in figure 5. 

The equations (1±) to (7) for stress, deflection, 
and acceleration give the values associated with the 
nth mode of vibration. Since all modes of vibration 
occur simultaneously, the net results are the superposition 
of the effects of all modes. This superposition gives 
the following equations: 

For bending stress, 



a(x,y,t) = - ^E\A x e 






(8) 



For average shear 



stress , 
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For deflection 

L 2( It 
w(x,t) = - p V^ie ^ sin li^t 

+ C 2 e 2E fc sin <o 2 t + •••/ ( 10 ) 
For acceleration, when damning is sufficiently small, 

a(x,t) = I f l^i^l^ s i n ^ 

+ C2^2 2e ^ S sin oopt + •'•/ 

The equation for bending stress (equation (i<)) 
reveals that, at a given percentage of the distance from 
root to tip, the bending stress for a particular mode is 
independent of the length of the beam and depends only 
on the velocity before impact. The equation for shear 
stress (equation (5)) reveals that the shear stresses at 
any station vary inversely with the length of the beam. 
These results are contrary to those that might be 
expected on the basis of experience with the static 
behavior of structures. For this reason an experimental 
investigation was made. 



Experimental 

A circular steel tube of 1-inch outside diameter 
and 0.028-inch wall thickness was mounted symmetrically 
on the end of a pendulum to form a pair of cantilever 
beams. (See fig. 6.) The pendulum was permitted to start 
its swing from a predetermined position and was suddenly 
brought to rest at the bottom of its swing against an 
electromagnet used to prevent rebound. The effect of 
length was studied by reducing the length of the tube in 
successive tests. The bending and shear strains were 
measured by electrical strain gages that were mounted on 
the tube as shown in figure 7. Each pair of gages was 
incorporated into a Wheatstone bridge circuit as shown 
diagrammatically in figure 8. The outputs of the bridge 
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systems were fed through a strain-gage amplifier into a 
multichannel oscillograph that recorded the strains on 
moving photographic paper. The amplitude of the components 
of strain due to the modes of higher frequency was 
reduced, however, because of the response characteristics 
of the oscillograph. The frequency-re snonse curve for 
the oscillograph used is given in figure 9. 

Typical records for tubes of two lengths are shown 
in figure 10. Inspection of the record for the cantilever 

beam 26^ inches long shows the superposition of the 

second and third modes upon the first mode. The record 
shows that, in the case of the bending strain, the 
contribution of the second mode is small; whereas, in 
the case of the shear strain, the contribution of the 
second mode is large. This observation confirms 
qualitatively the theoretical results shown in figure ij.. 
The same effect is not shown, however, in the record for 

the cantilever beam 11^ inches long because of the 

combined action of damping and reduced response of the 
oscillograph to the higher frequencies associated with 
this short length of tube. 

The bending stresses computed by use of equation (8), 
in which only the first three modes are used, are given 
by the solid-line curve of figure 11 for the cantilever 

beam 26jj* inches long. Comparison of this curve with the 
record obtained during the first \ cycle of the first 
mode (see fig. 10) shows good agreement as regards the 
wave shape . 

Because of the damping present in the tube and the 
response characteristics of the oscillograph, the only 
component of vibration that could be satisfactorily 
recorded for all lengths of cantilever tube was the 
fundamental or first mode. The quantitative results of 
the tests consequently were based upon this mode of 
vibration. This procedure is sound because the effects 
of the various harmonics are independent of one another. 
In the analysis of the results, the data had to be 
corrected for the influence of the magnet. 

The observed frequencies are compared with the 
frequencies computed from equation (1) for the first 
mode in the following table: 
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Frequency 


Length 

(in.) 


Observed 
( CDS ) 


Computed 
( cps ) 




17.5 


17.5 




27 .9 


28 .2 




52.1 


53.2 




131 


137 


_ "1 


272 


277 

1 



The experimental values of extreme-fiber bending 
stresses and the shear stresses at the root, for the° 
fundamental mode, are plotted in figure 12. in figure 12 
are also shown the corresponding theoretical curves of 
equation (1+) for bending and equation (5) for shear with 
n taken as 1. It is observed that the experimental 
points follow the trend of and lie close to the theoretical 
curves . 
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APPENDIX 
THEORETICAL DERIVATION 



General analysis . - Consider a beam of uniform cross 
section in equilibrium. If a portion of the beam is 
suddenly disturbed, as by a shock, in a direction 
perpendicular to its length, the beam is set into damped 
bending oscillations. The equation of motion for these 
bending oscillations is given by the differential 
equation (reference 1) 

EP 2 4 + x P 2^ ♦ I »f= = o (ad 



The damoing term Ap > is derlved on the 

assumption that the longitudinal damping force per unit 
area at any point on the cross section of the beam is 
proportional to the rate of change of longitudinal strain 
at that point. (See reference 2.) This type of force is 
analogous to ordinary viscous drag, in which the tangential 
force per unit area is proportional to the rate of change 

of shear strain. With the use of the notation c 2 = 
equation (Al) can be written 



64g k + 1 £ 2 w = 0 ^ A2 ^ 

6x^- E 6x^6 1 c 2 p 2 dt 2 

In accordance with the Heaviside operational methods 
(reference 3), equation (A2) may be reduced to an 
ordinary differential equation of the fourth order by 

ft 

writing p = — ; thus, 
6 1 
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The general solution of equation (A3) is 



W 



= P cosh 97 + Q sinh + R sin G7 + S co3 G~ 



L 



(A^) 



where 



9 = L 




The coefficients P, Q, R, and S are bo be determined from 
the boundary conditions. The case under consideration is 
that of a cantilever moving with uniform velocity v and 
having its base brought ins tanbaneously to rest. The 
boundary conditions for this case are 



The velocity of the root as given by the first boundary 
condition is represented graphically in figure 13(a). 
The rules of the Heaviside calculus, however, have been 
devised for a disturbance, called the unit function i , 
shown in figure 13(b). By the principle of superposition, 
the velocity function shown in figure 13(a) may be 
considered as a superposition of those shown in figures 13(c) 
and 13(d). The velocity therefore consists of a constant 
velocity v (fig. 13(c)) added to the solution of the 
problem obtained by the Heaviside expansion theorem for 
the disturbance shown in figure 13(d). On the basis of 
this procedure, the first boundary condition may be 
wri tten 
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With the application of the boundary conditions to 
equation (A4), the operational form of the solution for 
the velocity (that induced by the disturbance) is found 
to be 



4 



pw 2(1 + cosh 8 cos 9) 



(1 * cos 8 cosh 8) (cosh 9^ + cos 9~ ) 



+ sin 8 sinh 6 (cosh Qj - cos 



+ (sinh 9 cos 6 + cosh 9 sin 9) (sin 9^ - sinh 9^)J 



(A5) 



Interpretation of this operational expression and addition 
of the constant velocity v gives for the total velocity 



6w( x, t) 



vl + &r> F(9 nf7 )e 
ri=l 



cos OOn't 



2E sin V 



1 - 



^n 2 
iiE 2 




(A6) 



where 



'n 



nth positive root of 1 + cos 9 cosh 9 = 0 



8 2 



undamped natural angular frequency of 
nth mode, radians/sec 



co- 



il 



\ 2 C0 2 

n 



damped natural angular frequency of 
nth mode, radians/sec 



x sin 9 n sinh 9 n (cosh 6 n ^ - cos - (cosh 9 n sin 9 n + sinh R n cos 9 n ^ (sinh 9 n £- - sin 9 n ~} > 

V 1 !/ 0 n (cosh 9 n sin 9 n - sinh Q n cos 9 n ) 



w 



(*,t) - 2v ) ^**/o » sin » B 'ti 

I c J e""^ t sin u) n 'ti (A7) 



where 



2F ( e nr) 



e, 



n 

The contribution of the nth mode to the deflection is 



1 ^TT=r-t 

J- ^ H 



w n (x,t) =|^€ n r _ i_ e" sin % tti (A8) 



o 



as 
o 



Integration of equation (A6) with respect to the time with the condition (w) t=0 - 0 
gives for the deflection £- 

ro 

->3 



V>4 



11+ 
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When > 1, equation (A8) may be put in the form 

2E 



w n (x,t) 



_ v Lf 
c p 



\oo~.2 



'n 



v2(o 2 

*T1 



sinh o) n f t i ( A9) 



2 



where now 



|\2m 2 



oa n ' - u) n \ 



n 



I|E 2 



The form indicated by equation (A8), where — -~ < 1, is 

characteristic of the lower modes and represents damped 
oscillatory motion. The form indicated by equation (a9), 

\oc n 

where -5— > 1 (damping greater than critical), is 

characteristic of the higher modes and represents 
subsidence motion. 

From equation (Ao) for velocity and equation (A7) 
for deflection, the complete behavior of the cantilever 
may be determined. The quantities of interest are the 
bending stresses, the shear stresses, and to some extent 
the accelerations. When damping is present, the equations 
representing the contribution of the nth mode to these 
quantities may be given in the two forms indicated by 
equations (A8) and (A9) • In subsequent equations, however, 
only the form indicated by equation (AS) is given because 
it is characteristic of the modes that are of practical 
importance . 



Bending stresses The bending stresses o(x,y,t) 
at any fiber distance y from the neutral axis are 



A > 

o(x,y,t) - Ey ^— - 

dx* 2 fe 



fi | f ) A,. ] e 2E sin co n 'ti g- 



H 



where 

sin 9 n sinh 9 n (cosh 0 n ~ + cos ^n^) ~ (cosh B n sin 9 n + sinh 8 n cos 9 n j (sinh 9 n ^- + sin 9 n ^) 
A n = 2 - - 1 

9 n (cosh 9 n sin Q n - sinh 6 n cos 9 n J 



The bending stress due to only the nth mode is 



a n (x,y,t) = E * A n gg ; e 2E sin co n 'ti 

if - x ^ 



Ue 2 

H 



Shear stresses .- The average shear stress over the cross section 7(x,t) is 
T \ x , t ) - T 

6x3 



v — 



= E J £ > S n | 1 _ e 2E sin a) n 't i 



a., 2 



n=l I /V ^n 



l«s 2 

where 

sin Q r sinh Q n (sinh fi n ^» - sin 6 n r) - (cosh 9 r . sin Q n 4- sinh D R co s g n ^ (cosh G n ^- * cos G n ^) 
n cosh P n sin Qn - sinh 8 n cos 6n 

The average shear stress due to only the nth mode is 



f n (x,t) = E | £ Bn 1 e~ 2E ° sin o) n 'ti > 

J a - I 

t! 

p- — 

Accelerations • - From equation ( Ao) , with the aid of the relation m 

pP(t)i = F(0)pi + P T (t)i 
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the acceleration anywhere on the beam is found to be 



a(x ,t) = 6 2w ( x ,' t } = 



dt 



2 



n=l 



n=l 



p2 



e 2E sin 0) n ' t 



— V 

't + ] - 



UE2 



1 - 



2E 2 



cos a: n 't ]jl 



With the aid of the orthogonal properties of the 
functions x '( e n^) it is possible to show that the 

CO 

quantity 2 ^^ ? ( 9 n~) * 1 reduces to zero when 0 < - ^ 1 - 

n=l m 

At ^ = 0, the quantity 2^)T P vn^) ec l ua ls zero, and only the 



n=l 

term -vp-/ remains. This term indicates that at t = 0 
an infinite acceleration of zero duration exists at the 
root . 

The acceleration due to only the nth mode is 




if \ 1 - — 

f J4E 2 \ 2fi2 

Comparison with the expression for w n (x,t) (equation (a8)) 
shows that the acceleration for each mode is out of phase 
with the deflection. When damping is sufficiently small, 
however, the relation between the acceleration and the 
deflection reduces to the well-known result for undamped 
vibration 

a n (x,t) a -<*> n 2 w n (x,t) 
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Fig. 3 




Figure 3.- Variation of shear - stress 
coefficient B n with x/L. 
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Fig. 5 




Figure 5. - Variation or deflection 
coefficient C n with x/L . 
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Figure A?.- Graphic representation of 
various velocity functions. 



